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We investigate the time dependence of the quantum statistical model with A(^*/4! self-interaction 
and consider the resulting induced particle number density. For a special example in which the 
classical approximation exhibits parametric resonance, the effects of the back reaction and especially 
of collisions, treated in a perturbative way, are analysed. 

I. INTRODUCTION 

A widely used method for the investigation of the non-equilibrium evolution of quantum fields is the decomposition 
of the field into its mean value (or zero mode) and fluctuations, and the solution of the resulting coupled equations of 
motion. While the mean field rolls down some potential and/or oscillates around its minimum, its energy is transfered 
into the modes, a process resulting in particle production. This method has been applied to models of preheating in 
inflationary cosmology , to several examples of non-equilibrium dynamics of phase transitions ||^,|^ in the early 
' universe and heavy ion collisions, to the calculation of the formation of DCCs and also to more exotic processes 
like the decay |^ of parity odd metastable states in hot QCD ||T^ . 

In general, the coupled non-linear equations of motion which describe the real time evolution of the quantum fields 
cannot be solved exactly. The natural first step, then, is a classical approximation, in which the resulting mode 
equation exhibits parametric resonance for certain parameter values, leading to a huge amplification of particle 
production. This approximation is applicable for the early stages of many processes, when the number of particles 
produced is not too large. 

I In the next step, the back reaction of the created particles has to be included. This is usually done with a term 
that changes the effective mass in the time-dependent oscillator equation. Most existing calculations make use of the 
Hartree-Fock- (or large N-) approximation, in which the coupled equations for the mean field and the fluctuations are 
solved self-consistently. If and how fast the inclusion of the back reaction influences or even destroys the parametric 
resonance, depends strongly upon the model under investigation. For example, in certain models of preheating with 
small coupling, the back reaction is not very destructive whereas in the model for the decay of parity odd 

metastable states Q, particle production is strongly suppressed as compared to the classical case. Heuristically, 
this could be explained with the growing effective mass of the particles which makes them harder to be produced. 
Eventually, the Hartree-Fock approximation fails when the fluctuations grow comparable to the value of the mean 

; field. 

When the density of the particles created becomes large, the effects of collisions like scattering off or into the 
resonant modes, additional particle creation, dissipation to modes with higher energies and maybe thermalization, 
, are expected to become important. Parametric resonance amplifies only certain low momentum modes, producing 
highly non-thermal states with the energy concentrated in the infrared. In order to approach thermal equilibrium, 
energy must be scattered into higher momentum modes. These effects are not taken into account by the Hartree-Fock 
approximation, which corresponds to a truncation at the level of two-point functions. An understanding of these 
processes which are expected to terminate the parametric resonance and should eventually lead to equilibration are 
crucial not only for the scenario of inflationary cosmology (for a review see Q), but also for a deeper understanding 
of the different stages in heavy ion collisions 0. 

In the present work, we consider a field theory model with Jj^'* interaction and propose to include collision terms 
represented perturbatively by the sunset diagrams of O(A^), but expressed by the full two-point Green function in 
order to have an approximation applicable far from equilibrium. For the emergence of parametric resonance [0, 
typical parameters in the equation of motion must be periodically varying in time (e.g. the frequency in a harmonic 
oscillator equation), but the number of spatial dimensions is of no importance. Therefore, as a first step we consider the 
case of zero space dimensions, effectively a quantum statistical model, which simplifies the calculations considerably. 

After introducing the set of equations in section II, their solutions and the application to the preheating model of ||] 
are studied in section III. It is shown how the parametric instability of the classical approximation is destroyed when 
the collision term is switched on, i.e. when additional transfer of energy between the zero mode and the fluctuations 
is introduced. 
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II. QUANTUM STATISTICAL MODEL WITH COLLISION TERMS 



We consider a real one component scalar field with coupling, but for simplicity with zero spatial dimensions. 
Following the standard method as described, e.g., in |]3|-p|,p^, the field operator is decomposed into its zero mode 
ip = {4>{t)) and fluctuations x, 

m = ^{t)+m (1) 

with {x{i)) — 0. The full coupled, nonlinear equations of motion read 

m + m\{t) = -^v\t) - ^^{t) {x\t)) - A {^\t)) (2) 

for the zero mode and 

m + m^m = -^v^'mit) - ^ (x^ w - ix'm m - ^ (x^w - ix'm (3) 

for the fluctuations, which is an operator equation in x- In the following, we will make use of the mode functions 
x(i), which are introduced in the same way as in the usual Fourier decomposition of the Heisenberg operators in field 
theory: 

i{t) = ^{ait)+a\t)) (4) 
^ {ax{t)+a^X*{t)) (5) 



'0 



where ujq = w(to) is defined below in (15). The time- independent creation and annihilation operators a = a(to), a* = 
a^{to) act on the initial Fock vacuum state. The density matrix at to is chosen to fulfill (a) = (a^) = 0, {aa) = 
(a^a^) = etc., except for (aa^) and (a^a) which are assumed to be nonzero. 

The first step towards the solution of the equations of motion is the classical approximation, which is linear in the 
mode functions and does not include any back reaction effects of the particles created: 

^ + mV(t) + ^v'(t) =0 (6) 

X+(ni' + ^^^t)^x{t)^0- (7) 

(For the solution, see section III.) In this case, the induced particle number density, which is defined as 

nit) = {aHt)ait)) (8) 
_ Trat(t)Q(t)pfa) 

Trp(to) ^ ' 

^ TVat(to)afa)p(t) 
Trp(io) 

can be written in terms of the mode functions as 



(10) 



r^W = 7 lxWP + ^ (11) 



1 , \m\'\ 1 

2 

just as for the harmonic oscillator, with the initial conditions x(^o) = 1 and x(^o) = — *'^(io) consistent with n(to) = 0. 

A consistent method to include back reaction effects of the fluctuations on the mean field and on themselves is the 
Hartree-Fock approximation x'^ 3(x^)Xj leading to (x^) = 0. The equations of motion are 

m + m\{t) + A^3(i) = _^ (^^^) ^{t) (12) 

m + m\{t) + ^^\t)x{t) = -\ ix') X(t) (13) 
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with the back reaction term 

{X^)^7^\x{t)\\ (14) 
In this approximation, the Hamiltonian is stiU that of an osciUator, but with the time dependent frequency 

u.'{t)=.i'+\v\t) + \{ie). (15) 

The definition of the particle number in this approximation can be kept as before (for a discussion of the time- 
dependent oscillator see, e.g., [^), and the order in A is the same. 

Although the mode functions will be used below for the numerical calculations, we also need their relation to the 
Green functions to establish the equations of motion with collision terms. The Green functions are defined as (see, 
e.g., 0) 

G>M = {mx{t')) (16) 
and can be written in terms of the mode functions as 

G>{t,t') = [(1 + n(io)) X{t)x*{t') + n{t^)x*{t)x{t')] . (17) 

ZUJQ 



With n(to) = 0, this simplifies to 



G>{t,t')^^x{t)x*(t'), (18) 



and G^{t,t') — G-^{t',t). Because we use the mode functions for further calculations, we effectively only need the 
equation of motion for the equal-time Green functions 



[d^ + m^ + j:\t)]G<{t,t%,=t^~z I dt"J:>{t,t")G<{t",t%,=t (19) 

' dt"j:<it,t")G>it",t%,^u 



which is the Schwinger-Dyson equation (see, e.g., ||lj|) in the limit of equal times t' = t (for the notation, see IfiSl). 
The integration contour is shown in fig.l. 
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FIG. 1. The Schwinger-Keldysh closed-time-path contour for the case t' — t. 
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The tadpole self-energy 



(20) 



is the part which already appeared in the Hartree-Fock approximation ([12[),(|13|). In order to close the system of 



equations by including effects due to collisions, we approximate the self-energy in eq. ( 
function (|l8|). In the presence of the condensate field (p{t), the sunset-type approximation (cf 



19| ) by the two-point 



is generalized to (see fig. 2) 



(21) 



S>(t,i") 



A2 



[G>{t,t")]' - \^it) [G>{t,t")\\{t") - ^^^(t) [G>{t,t")\ ^^{t"). 



(22) 




FIG. 2. Graphical representation of the r.h.s of the Schwinger-Dyson equation for the Green function G(t,t') (solid curve). 
The dashed curve denotes the mean field. The statistical factors (cf. [16]) are given. 



Expressing the Green functions in terms of the mode functions, we obtain the equation of motion up to O(A^): 



x(i)+(m' + ^^'(i) + ^lxWnx(0 



A2 



24^o' J to 



dt'Im 



H J to 



4cj, 



dt'lui 



dt'lm 



Jto 



^{mt')(x{t)x*{t')r 



(23) 



where 



(24) 



Correspondingly, the equation of motion for the mean field is (see fig. 3) 



d^,+m' + ^^^{t) + ^G{t,t) 



V{t) = -i / dt'{t>{t,t')-t<{t,t'))^{t') 



iX^ / dt'(p{t') 



{[G>{t,t')f-[G<it,t')f^ 



+-^it){[G>{t,t')f-[G<{t,t')f) ^{t') 
+^v\t) {G>{t,t')-G<{t,t'))^^{t') 



(25) 
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and the initial conditions are chosen as f(to) = ipo and (^(to) = 0. Inserting the mode functions, eq. (|2q) reduces to 



where it should be noted that ip{t) is a real function. 



A2 



24^o' Jt, 



dt'lm 



v(t') ixit)x*it')f 



(26) 



Jto 

A2 



dt'lm 



dt'lm [vp'(t)^^(t') {X{t)x*it'))] 




FIG. 3. Graphical representation of the r.h.s of the Schwinger-Dyson equation for the mean field (dashed curve). 



The particle number including collisions is calculated in general from 

h{t) =i[H,n{t)]. (27) 
Following the derivation given in ||l^,|l8| , this reads for a scalar field for the contributions due to collisions 

(28) 



where t and t' are on the contour of fig.l. Keeping terms of O(A^), we write 



^0 Jto 



d 



{t)coii^' — [l dt"^>{tX)-^G<{t",t')+ I dt"E<(t,t")^G>(t",f 



to a 



dt' 



(29) 



where S is given by the graphs of fig. 2. With G expressed in terms of the mode functions, the result is, after integrating 
and adding the contribution ( |ll| ) without collisions, 



(30) 



dt' 



dt"Im 



to 



8^0 Jto 



^ / dt' 

8^o' Jto 



dt' / dt"lm 

to 

't' 



^{t'Mt")x{t')x\t') {x*{t")f 



dt"lra 



to 



^'{t')^'{t")x{t')x{t'){x*{nr 



(31) 



In the case of vanishing condensate field, i.e. ip — 0, but with collisions, the particle number in leading order of A is 

A2 sm^{2uja{t-to)) 



ncoii{t) 



48^4 



8^2 



(32) 



which is easy to find after approximating the mode functions in terms of free fields, x(^) = exp(— iwoO- This 
corresponds to the expression (3.21) in for the d = case, neglecting 0{X^) contributions. Here, we discuss n{t) 
of eq.(pO|) without expanding the mode functions in terms of the free field solution, but instead determine x(t) from 
the self-consistent set of coupled equations (ESh and (p§). 
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III. RESULTS FOR MODELS WITH COLLISION TERM 



In order to study the effects of the collision term, we solve the coupled equations of motion numerically for the 
three steps of approximation. 

As a first example, we consider the set of parameters — 1.6, A = 0.1 and (po — 7.75, and the other initial 
values for (p and x ^-re as given as in section II. (This choice will be justified below.) The results for the classical 
approximation eqs.(|^) and (0) (solid curve), including back reaction (dashed curve) and additionally including the 
collision term (dot-dashed curve), are shown in fig. 4. The growth of the induced particle number density does not 
exhibit parametric resonance for this set of parameters, but it is sufficient to illustrate the dramatic effect of the 
collision term: The inclusion of the back reaction (eqs.(p^,(|T^)) already suppresses this growth, and the additional 
inclusion of the collision term (eqs.(|2^),(p6|)) destroys it completely. 

3.5 I , 1 1 1 1 
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FIG. 4. Particle number density for = 1.6, A = 0.1; for the classical approximation (solid curve), including back reaction 
(dashed curve), and additionally including the collision term (dot-dashed curve). 
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FIG. 5. Zero mode for = 1.6, A = 0.1; for the three cases as in fig. 4. 
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As the main example, we apply our approximation to a special preheating model of chaotic inflation treated in 
where the development of parametric resonance does not depend on the expansion of the universe. The equation for 
the zero mode 

¥>"(a;)+^3(a;) =0, (33) 

written for the conformal field if — acj) (a is the cosmological scale factor) and the dimensionless time x — ^/Xiporj {rj 
being the conformal time) has the solution 



ip{x) =ipQcn[ x, 



V2j ' 

a Jacobian cosine function, leading to a Lame equation for the fluctuations of ip, 



.2 + 3''^^^)^ 



(/3k (x) = 0, 



(34) 



(35) 



where the comoving momentum k appears in = k'^/XipQ. The Lame equation is known to have unstable solutions 
of the form ip^ix) = P{x) exp(/i„a;), where P{x) is a periodic function, for certain values of k which form the so called 
instability bands. The Floquet exponent characterizes the growth of the solutions and of the particle number 
density In n„ « 2fi^x, and is real for the instability bands. In this special case, the instability band is given by 
3/2 < < V3 with its maximum value fJ,K.,max — 0.03598 at k^^^, ~ 1.615 As discussed in ||l^, the inclusion of 
different "back reaction" /rescattering effects will change the equations of motion, the frequency of oscillation and the 
effective masses. The instability bands are restructured, and scattering may lead to additional particle production 
and removal from the resonance. 

To investigate effects of this kind with our approximation, we first have to introduce the dimensionless time variable 
T = y/X/6Lpo t. Comparing our set of eqs. (|^), (0) with the corresponding equations of leads to the condition 



6m 



(36) 



i.e. the momentum k, which does not appear in our equations because we consider zero spatial dimensions, is 
"replaced" by the mass. Moreover, we have to neglect the mass term in the zero mode equations (|^) , (|l2|) , and (|26|). 

We choose the initial value ipQ — •y/6/A, which has the advantage that t ~ t, i.e. our time scale is the same as in 
1^. Choosing A = 0.1 and considering the special value n'^ax — 1-615 fixes m? for the mode equations. (This set of 
parameters corresponds to those of the first example.) Fig. 6 shows the results for the particle number density, where 
the classical result (solid curve) reproduces that of fig. 3 in The solid curve in fig. 7 corresponds to the solution of 
eq.(|3^). The parametric amplification of the classical approximation is suppressed only weakly by the back reaction. 
In contrast, the resonance is completely destroyed by the collision term given in eqs.(|23|) and (26). 




FIG. 6. Particle number density for m? = 0, A = 0.1 and — 1.615; for the classical approximation (solid curve), including 
back reaction (dashed curve), and additionally including the collision term (dot-dashed curve). 
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FIG. 7. Zero mode for = 0, A = 0.1 and = 1.615; for the three cases as in fig. 6. 



IV. CONCLUSION 

The essential feature of this model for the inclusion of collisions is that the two-point functions are not expressed in 
terms of planar waves, i.e. is not perturbatively expanded in terms of A. This is crucial for obtaining parametric 
resonance, which leads to amplified particle production, in the classical equations and even when the back reaction 
is included. With the inclusion of collisions, the parametric amplification disappears, at least within our treatment. 
But, unfortunately, our approximation breaks down already for small times, </ ~ 11 in this case (see figs. 6 and 7). 
Decreasing the coupling A, without leaving the resonance band of the classical solution, only shifts the instability of 

the approximation scheme towards larger times if — \J X/\tf. There is strong evidence that the effective potential 
becomes concave for t > tf, as can be seen from the evolution of the zero mode (fig. 7), which does not only grow larger 
than its initial value, but even seems to tend to infinity. Contrary to studies of thermalization (see, e.g., [po| ), we are 
not able to see the long-time behaviour, although the model includes (memory) terms which are non-local in time. A 
definite statement about equilibration, i.e. about particle number densities including momentum dependence, would 
additionally require the inclusion of d > spatial dimensions. 
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